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The quantum fluctuation of the relative location of two (n — l)-dimensional de Sitter branes 
(i.e., of n spacetime dimensions) embedded in the (n + l)-dimensional anti-de Sitter bulk, 
which we shall call the quantum radion, is investigated at the linear perturbation level. The 
quantization of the radion is done by deriving the effective action of the radion. Assuming the 
positive tension brane is our universe, the effect of the quantum radion is evaluated by using the 
effective Einstein equations on the brane in which the radion contributes to the effective energy 
momentum tensor at the linear order of the radion amplitude. Specifically, the rms effective 
energy density arising from the quantum radion is compared with the background energy 
density. It is found out that this ratio remains small for reasonable values of the parameters of 
the model even without introducing a stabilizing mechanism for radion, although the radion 
itself has a negative mass squared and is unstable. The reason behind this phenomenon is also 
discussed. 



Based on the idea of a brane- world suggested from string theory [jj], Randall and Sundrum proposed an 
interesting scenario that we may live on either of two boundary 3-branes with positive and negative tensions 
in the 5-dimensional anti-de Sitter space (AdS) [§J||. One of the attractive features of the Randall-Sundrum 
(RS) scenario is that the gravity on the brane is confined within a short distance from the brane even for an 
infinitely large extra-dimension . This applies to the positive tension brane and it is because the AdS bulk 
on both sides of the positive tension brane shrinks exponentially as one goes away from the brane. 

Since the RS scenario gives an exciting, new picture of our universe, it is clearly important to study various 
aspects of this scenario to test it or to give constraints on its parameters. One good example is an analysis 
done by Garriga and Tanaka 0, in which they have shown that the radion in the original RS scenario acts 
as a Brans-Dicke scalar on the branes at the linear perturbation order and the effective gravity is that of a 
Brans-Dicke theory with positive and negative Brans-Dicke parameters on the positive and negative tension 
branes, respectively, where the values of the Brans-Dicke parameters are determined by the distance between 
the two branes. In a previous paper, we have shown that essentially the same situation arises in the case of two 
de Sitter (dS) branes embedded in the AdS bulk. However, we have also shown that the radion effectively has a 
negative mass squared with its absolute value proportional to the curvature of the dS brane, hence is unstable 
if it can fluctuate by itself without the matter energy momentum tensor. 

The phenomenon we shall study in this paper is the quantum fluctuation of this mode, called the quantum 
radion. To make clear what we mean by the quantum radion, let us describe the mode in more detail. Our 
brane universe may be displaced from the Oth order trajectory of a homogeneous and isotropic brane. By 
appropriately fixing the coordinate gauge, the displacement perpendicular to the brane can be described by a 
scalar function on the brane. Assuming there are two branes that are fixed points of the Z<i symmetry, it can 
be shown that only the relative displacement of the two branes is physical, which we call the radion. There 
are two distinctively different kinds of displacement of a brane: the "bend" and "fluctuation" . The bend is a 
type of displacement due to inhomogeneitics of the matter energy-momentum tensor on the brane. The trace 



I. INTRODUCTION 



of the energy-momentum acts as an additional tension of the brane, and the brane must "bend" accordingly. 
The relative bend is described by the mode of radion that couples with the source on the branes, and it can 
be written as a functional of the energy-momentum tensor . In the context of the quasi-localized gravity 
discussed by Gregory, Rubakov & Sibiryakov S , the role of this type of radion has been extensively studied 
[iTHnj. On the other hand, the "fluctuation" is a type of displacement that is purely geometrical, which obeys 
a free wave equation without source. The relative displacement of this kind is the mode of radion which we 
shall discuss in this paper. This mode of radion was studied first by Charmousis, Gregory and Rubakov for 
the RS branes whose effective radion mass is zero by solving the field equations for the RS branes [fl2f and by 
Chacko and Fox Jl3| for the dS and AdS branes whose effective radion mass squared are negative and positive, 
respectively. 

The fact that the radion mass squared is negative (or zero) suggests the (marginal) instability of the two 
brane system. In fact, in the case of the original RS flat two brane model ||, the negativity of the Brans-Dicke 
parameter on the negative tension brane Q] can be regarded as a result of this marginal instability. To recover 
the stability, Goldberger & Wise introduced a bulk scalar field (lj] that couples to the branes in such a way 
that the distance between the two branes is stabilized. However, it should be also noted that the effective 
Brans-Dicke parameter on the positive tension brane is positive and it can be large enough to be consistent 
with experiments for the separation of the branes larger than the AdS curvature radius at least at the linear 
perturbation order [||. Hence a stabilization mechanism may be unnecessary if we live on the positive tension 
brane. We therefore do not introduce a stabilization mechanism. 

We shall work on the system that consists of an (n + I)-dimensional AdS bulk spacetime bounded by two 
branes of constant curvature that are fixed points of the .^-symmetry. The zero-curvature branes correspond 
to the flat RS branes [^|J3| , while the positive-constant curvature branes correspond to the dS branes |l5| . Our 
main concern is of course the dS brane case, but we treat the flat brane case simultaneously to make clear the 
similarities and differences between the two cases. 

The dS brane case is of particular interest because it gives a good model of braneworld inflation. In the 
standard 4-dimensional inflation, the quantum vacuum fluctuations play a very important role. It is therefore 
natural to ask if the quantum radion fluctuations play an important role, if not disastrous, in the braneworld 
inflation. 

It should be mentioned that historically a very similar situation was analyzed by Garriga and Vilenkin 1 16 1^] 
in which they considered the fluctuations of a thin domain wall in (N + 1) spacetime dimensions. Although 
they assumed the Minkowski background, many of the results obtained there apply equally to the present 
case. In particular, they showed that the wall fluctuation mode is represented by a scalar field living on the 
A-dimensional de Sitter space which describes the internal metric on the domain wall, and the scalar field has 
the negative mass squared —NH 2 . 

The paper is organized as follows. In Sec. [fl], we describe the background spacetime and our notation. In 
Sec. [ill) , we solve the perturbation equation in the bulk that describes the radion mode. We find there is a 
gauge degree of freedom that should be carefully treated in the dS brane case in contrast to the flat brane case 
where no such subtlety arises. In Sec. |Fv|, assuming the dependence on the extra dimensional coordinate that 
solves the perturbation equation in the bulk, we derive the effective action for the radion and quantize it. In 
Sec. |v|, based on the result obtained in Sec. IV, we evaluate the effective energy density of the radion on the 
brane which is present at the linear order in the radion amplitude, and estimate its effect by calculating the rms 
value. We find the effect remains small for reasonable values of the model parameters even though the radion 
itself is unstable. In Sec. VI, we summarize our results and discuss the implications. 



II. BACKGROUND 



First, we summarize the basic equations of the system and explain our notation. The system we consider is 
a Zi symmetric (n + l)-dimensional bulk spacetime with two (n — l)-dimensional branes as the fixed points of 
the symmetry. The bulk metric g a b obeys the (n + l)-dimcnsional Einstein equation: 



G a b + h- n +l9ab — K 2 T a b , (2.1) 



where n 2 is the (n+l)-dimensional gravitational constant, A„ + i is the (n+l)-dimensional cosmological constant 
which we assume to be negative, and T a b is localized on the branes. We use the Latin indices for (n + 1)- 
dimensional tensor fields in the bulk and the Greek indices for n-dimensional tensor fields on the brane. Denoting 
the rt-dimensional metric on the brane by q^ Ul we decompo se th e localized energy momentum tensor T a b into 



the tension part —aq^ and the matter part t^. Then Eq. fl2.l| ) reduces to 
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{n) G^ + A n q^ = K 2 n T„ v + K 4 ir^ - Ep,, (2.2) 

on the brane |[18| , where 7r^ w is a tensor field quadratic in t^, and is the projected (n + l)-dimensional Weyl 
tensor defined by :— ( n+1 J 1 C a ^ lbu n a n b with n a being the unit vector normal to the brane. The constants A„ 
and are related to the basic constants of the systems asf] 

A„ := ^ (A, i+1 + «V) , := ^-^ (2.3) 

It may be noted that although the decomposition of T a (, is not unique, there is no arbitrariness in the metric 
on the brane. 



We consider branes with = 0. Looking at Eq. (2.2), we see that the E^ term can be considered as the 



effective energy-momentum tensor on the brane induced by the bulk gravitational field, 

<:=-i^- ( 2 - 4 ) 

The background we consider is the RS brane system and the dS brane system. To be more precise, the bulk is 
the (n + l)-dimensional anti-de Sitter spacetime (AdS n +i) whose metric may be written as 



ds 2 = g ab dx a dx b = b 2 (z) [dz 2 + ^ v dx fi dx u ] ; b(z) = : ^ , I =- " ' " ' ' 

sinh V K-z 



2A 



n+l 



(2.5) 



where 7^ is the metric of a Lorentzian n— dimensional constant curvature space with curvature K, = or K, = 1, 
and t is the curvature radius of AdS n _|_i. The function b(z) is called the warp factor. The two branes are placed 
at the coordinates z = z + and z = z_ (z + < z_) with their tensions and cosmological constants given by 

,2(n-l) , . r- . , , 1 (n- l)(n- 2)/C . . 

<r± = ± \ 2£ ; cosh(V^z ± ) and A n , ± = ^ ; , (2.6) 

which satisfy the effective Einstein equations with T^„ i= t = an d E^± — 0. The choice K, — corresponds to 
the RS brane system while K, = 1 corresponds to the dS brane system. Note that the n-dimensional gravitational 
constant k„ on the positive tension dS brane is related to that on the positive tension RS brane, say rs , as 

K n = K ^RS COSll2; +- 

Although we perform most of our calculation in the coordinate system with z, sometimes it is more convenient 
to use the coordinate system with the proper distance coordinate r defined by 

dr = -b(z)dz. (2.7) 

It is also useful to introduce the rate of the change in the warp factor b{z) in the r coordinate, which we shall 
denote by J(z), 



d r b d z b coshv^z 



III. RADION MODE IN THE BULK 

In this section, we consider the gravitational perturbation on the background described in the previous section. 
Our analysis is a generalization and a reformulation of p2l which discussed the RS brane system (fC = 0). The 
4-dimensional dS and AdS brane system was discussed in JT3J . 

We denote the gravitational perturbation by h a b to the background metric g a b, i.e., g a b = gab + h ao where g a b 
is the full metric. We choose the RS gauge with respect to the positive or negative tension branes and denote 



'There is a typo in the corresponding equation in 
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the gravitational perturbation in this gauge by ftj^j and the coordinates by {T[±] , ^u-i}- The RS gauge condition 
is 



In the RS gauge, the Einstein equations (2.1) simplify to 



—bd r b n ~ l bd r + U r 



IK 



h [±] 



(3.1) 



(3.2) 



Since h^n a = 0, we may use Eq. (A23) of (|] to obtain the boundary condition on /iS . Expressing the 
perturbed locations of the branes in terms of ip^ (x) as 



n+i = r+ + ip^(x) and n_i = r_ + ip^ '(x) , 
where r± denote the positions of the background branes, we have 



h [±] 
If 



at r [±] = r± =F , 



where 



(3.3) 



(3.4) 



(3.5) 



Note that t^± is set to zero as we are interested in the gravitational perturbation that does not couple with 
the matter on the brane. It should be noted here that Eq. (A23) of || is defined from the side with the smaller 
value of r, while in this paper the boundary condition at the negative tension brane is g iven from the larger 
value of r; rr_i = r_ + 0, hence the signature is reversed. In passing, we also note that Eq. (3.4) may be derived 
by first introducing the Gaussian normal coordinates with respect to the branes and transforming them to the 
RS gauge, as was done in [Ejfl. 



With the junction condition (3.4) at hand, we consider a particular solution of the form, 

The traceless condition on the metric perturbation in the RS gauge demands tpktf to satisfy 

[-□„-n/C]^±] =0, 



which implies 



[-□„+n/C]£^ [±] =0. 



Then Eq. (gj) reduces to 



—d z b n - L d z + {n- 2)K 



u(z) = 0. 



(3.6) 
(3.7) 
(3.8) 

(3.9) 



The general solution to this equation is obtained as follows. First, it is easy to see that u^(z) = U(z) is a 
solution. Note here the identity d z J(z) — K.b^ 1 (z). Then the other independent solution, u^(z), may be found 
using the Wronskian of the above equation, which is 



W(z)=£ n - 2 J(z) / dz' J- 2 (z')6 1 - n ( 2 ')' 



Thus, the general solution of the form (3.6) is 
where and flW are constants. 



(3.10) 



(3.11) 
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Note that there is an ambiguity in u^(z) which depends on the choice of the integration constant. The 
change of it induces the change of the coefficient C™. This implies the solution represents a gauge degree 
of freedom. In fact, if we evaluate E^, , we obtain 



(3.12) 



which is independent of . Note also that the z-dependence of E^ is clearly consistent with the boundary 
condition at the branes with no source on it, d z (b n ~ 2 E^ u ) = ||. 

Since E^ v is a gauge-invariant tensor field, we see that the displacement of the branes should be related as 



(3.13) 



With the above relation at hand, we can use a residual gauge degree of freedom in the RS gauge to bring the 
two different sets of the coordinates {a^j^j} to the same coordinates, say {x a }. This is done by an infinitesimal 
coordinate transformation, 



r-»r + f(^), 



c ^ + is- 2 f^(/) + f(/) 



£r = % 

S 2 



CM - 



[±] 



where a is the yet remaining gauge parameter and the function B{z) is defined as 

f b- 2 (K = 0) 



b- z dr 



2 cosh V/Cz - 1 



I 2 K 

The metric perturbation in the unified coordinates {x a } is 

V = bH~ l (au^(z) + u< a >(*)) flW^W , 
and the branes are now placed at 



(J -I- 1 ) (AC = 1) 



r± + 1 + — 



[±] 



(3.14) 
(3.15) 

(3.16) 

(3.17) 
(3.18) 



We note t hat the positions of the two branes cannot fluctuate independently in this unified gauge because of 
Eq. (3.13). This shows the nature of the radion that it describes the relative displacement of the branes. 

With the above form of the metric perturbation and the locations of the branes, we find that the boundary 
condition on h a b gives 



£>[±] 



2f 



alC£ + b±d z u(. 2) (z ± ) 



1 + 7 



which implies 



1 + 



K 



2/' 



(3.19) 



(3.20) 



An apparent complication in the case of K. = 1 is that this seems to constrain the values of which is just 
unphysi cal, gauge parameters. This can be resolved by considering the integration constant term in vS 2 > defined 
in Eq. ( pMC| ). For a change of C^ as 



we consider the simultaneous change of the integration constant term in as 

_ u (2) _p[±]/3) u (i) , 



(3.21) 
(3.22) 
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which results in 







(3.23) 



This combined change of and the i ntegr ation constant term in u^(z) leaves Eq. ( 3.20 ) as well as h]^} 
unchanged. One can also see from Eq. (3.23) the reason for the simplification in the case of K = 0, i.e., the 
change of does not affect d z w- 2 ' if K = 0. 
Introducing (p{x) defined as 



£>[+y +1 



21 



27" 



the results are summarized as 



(3.24) 



• h 



flis 



2& 2 (cra (1) +uW)C ltv $ and = 2(b/lf- n C llv (p 



,(2) 



positive 
negative 

(2)f 



brane at r = r± + L(z±; a)(p 



where L(z; a) :— alK, + b{z)d z vS > {z) and <p satisfies 



nK]<p{x») = 0. 



(3.25) 



(3.26) 



As it should be clear by now, the parameter a describes the residual freedom of the RS gauge in the unified 
coordinates, corresponding to the coordinate transformation, r — ► r + /C/(x M ) and x^ — > x^ + Jf^(x^) with 
f(x») satisfying £^f(x») = 0. 



IV. QUANTUM RADION 



Gi ven t he initial data for tp, we obtain the full time evolution of the fluctuations of the brane by solving 
Eq. (3.26). As the initial data, it is natural to assume that <p> is classically null, and only the quantum vacuum 
fluctuations are present. 

To consider the quantum fluctuations of the brane, we need the action for ip in the first place. We obtain 
it by substituting t he cla ssical solution (3.25) into the original action of the system, but without constraining 



<p(x) to satisfy Eq. (3.26). This is because thus obtained action is already maximized except for the degree of 
freedom for the brane to fluctuate and hence is adequate for the quantization of such a degree of freedom . 

Our system consists of identical bulk spacetimes Mi and M2 with branes £+ and £_ that are the boundary 
hypersurfaces of Mi and M2 with tensions of <r+ and cr_, respectively; dMi = £+ U £_ and DM2 = £+ U 
The action of such a system is given by 



Ir-2K,x '■= r-j / V^gdrd n x 

J Mi 



I — Ir-2A,1 + Ir-2A,2 + Ir,1 
_ 2 A, 



I 



K.2 



n+l 



+ 1*^ 



Ia- 
1 



-=■ / ^f~qd n xK, I 



a Mt 



-qd n xa±, (4.1) 



where q a b '■— 9ab — n>a,n-b, K a b := q a c q b V c n>d and n a is a unit vector normal to th e branes, pointin g fro m Mi to 
M.2- The substitution of and the locations of the branes given by Eq. (3.25) into the action (4.1), We find 



& [0[ = 



(n- l)(n-2)£ r > 



where 



N, 



(n-2)K. •= 



071-3 



dz'b 



'1,3— « 



(n - 2)1 



d^x^f^i (p [□„ + nK\ (p , 



b-d z u( 2 \z-)-b + d z uW(z+) 



(4.2) 



(4.3) 



The normalization factor Nt n _2)K. IS found to be identical to the one introduced in || for the normalization of 
the radion mode of -E M „. Details of the derivation and its comparison with those by Chiba [^0| and by Chacko 
and Fox |li| are given in Appendix [X| 
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Since A r („_ 2 )AC > 0, we may normalize (p as 
^(s") := 



2{n-l){n-2)N {n _ 2)K l ( I 



r- <PW), 



(4.4) 



such that the action reduces to that of a scalar field ip{x) with the mass squared m 2 := ~nJCb + 2 on the positive 
tension brane: 



q+ d^tpdui/} 



2 / 2 

-m ip 



(4.5) 



where g+^i, := b^j^ is the intrinsic metric of the positive tension brane. Here and below, we focus on the 
positive tension brane, but the parallel arguments apply to the negative tension brane except for the difference 
in the normalization factor. 

In terms of ^(x^), we see from Eq. ( $.25 ) that E^ v is given by 



E„ 



(4.6) 



where 



N E 



-TL 



V (n-l)(n-2)N (n _ m i 

r 7/-^ -.a/3 r 



lt>,v 



□, 



(4.7) 



Here it may be appropriate to make a comment on the form of Eq. ( }4.6[ ) in which the operator C^ v is replaced by 
its traceless part. In the previous section, we considered only the classical solution for the radion fluctuations, 
hence we were allowed to use the equation of motion to express in terms of the radion field, and the 
equation of motion for the radion was derived from the traceless condition of the RS gauge. When we consider 
the quantum radion, however, E^ must be traceless at the off shell level as well. This means that the traceless 
condition of the RS gauge should not be used to express E^ v in terms of the radion field. Instead, one has to 
go back to the original definition of the Weyl tensor in terms of the metric per turb ation whose trace may not 
be null. This is the reason for the appearance of the traceless part of C in Eq. (4.6). 

Now that we have obtained the action for the fluctuation of the brane, which is of the form of the scalar field 
action, we may canonically quantize it in the ordinary way. Adopting the flat coordinate chart on the brane, 
which is specified as 



dsi 



b\ lllv dx^dx v = a 2 (r]) [-dry 2 + {dxf 



a(r]) 



(1 - K) - ICrj 



we obtain 



#7,cc) = J d n l k [a k i> k (r))e ik - x +a{^t{ri)t 



— ik-: 



where a& and d^. satisfy the commutation relations, 



[fife, o fe/ ] = S(k - k') , [fife, fife/] = [fi£, a\.,] = , 



and {ipk(v)} satisfies the equation of motion, 

d 2 _ (n-tyC d 
drf 



n 



df) 



, 9 nJC 
k ~ — 



with the normalization condition, 



2i 

,n-2 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



(4.12) 
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As a natural choice for the vacuum we choose the Bunch-Davis vacuum, for which we have 



-ikr/ 



{K = 0), 

^&; ( "-n-^) ( "- 1)/2 <l «*(-*»>) (ic = i), 



(4.13) 



where fl^ 1 (z) is the Hankel function of the first kind. For K, = 1, this would correspond to choosing the de Sitter 
invariant vacuum if the mass squared were positive. Note that the rms amplitude of the vacuum fluctuations 
per logarithmic interval of k exhibits an infrared divergence in the limit rj — > — 0, 



1 



\krj\ 



(4.14) 



We now focus on the case JC = 1, which is of cosmological interest, and consider the effect of the quantum 
radion on the brane. As noted, ip shows an in stability which grows as (— 77) -1 oc a(rj). However, what on e ca n 
detect on the brane is E^ v given by Eq. (4J3) which acts as the effective energy- momentum tensor, Eq. (2.4). 
To see its behavior, let us spell out the components of explicitly. The components for each k mode are 



£00= N E 



E 0j = N E (-k) 



1 



d_ 

77 \drj 

— + - 

drj 77 



1 

i] 



Ea= Nsk^kYij + N E 6i 



-(- 



where Eq. (4.11 ) has been used, and Y, Yj, and Yij stand for 



Y — e ife x 



Y < = -H Y - 



n - 1 



Y . 



(4.15) 



(4.16) 



Inserting the explicit form of ipk given in Eq. (4.13), we find 



-^00 = 


N E k 2 {- 


-kr])( n - 


-l)/2o-(l) ,_ 
/3 (n-3)/2V 


-krf) Y , E% 


E j = 


N E k 2 {- 


-kr])( n - 


-l)/2rr(l) /_ 
a (n-l)/2\ 


■kr))Yj, 


E ij = 


N E k 2 {- 


-krj)^ 1 - 


-l)/2w(l) (_ 


-kv)Yij, 



n-1' 



-E, 



00 ■ 



where E^ and EE denote the trace and traceless part of Eij , respectively, and 



N E = N, 



(n-2)/2 



2 fc(n-l)/2 



(4.17) 



(4.18) 



From the above, and noting H^\z) ~ z~ u for z — > 0, we see that E^, which represents the effective anisotropic 
stress, dominates over the other components for k\rj\ — » 0, i.e., on super Horizon scales, though its physical 
amplitude (Ef^/a 2 ) decays as 1/a. 

As for the effective energy density induced by the radion, it decays rapidly as a -3 . However, it may be said 
that the decay rate is not fast enough in the following sense. Since E^ is traceless and conserved at the linear 
perturbation order, one can regard it as the energy momentum tensor of a radiation-like fluid. If it should 
behave exactly as radiation, the energy density would decay as a~ n instead of a -3 . Of course the decay rate 
a~ 3 does not mean E^ behaves as a dust matter. The fact that the decay rate is slower than the standard 
radiation fluid (for n > 4) is caused by a strong anisotropic stress present in the case of radion. Thus, in this 
sense, the instability of radion indeed shows up on the brane, but the effect turns out to be rather mild because 
of the nature of de Sitter space in which the background energy density remains constant in time. We defer the 
detailed discussions on this point to Appendix |c| 

To summarize, we conclude that the instability of the quantum radion induces a large anisotropic stress on 
the brane, but because of the rapid expansion of the brane it is not strong enough to give rise to a gravitational 
instability on the brane. 
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V. COSMOLOGICAL EFFECT OF QUANTUM RADION 



In the previous section, we have seen that the quantum radion fluctuations do not cause a gravitational 
instability on the brane. However, this does not necessarily mean the effect is always small. In particular, if 
the rms value of the effective energy density due to the radion field is comparable to the background vacuum 
energy density of the de Sitter brane, the evolution of the brane will be significantly affected. If we regard 
our background to represent a braneworld undergoing cosmological inflation, this would imply a significant 
modification of the braneworld inflationary scenario. In this section, assuming the positive tension brane 
describes our universe, we consider the cosmological constraints arising from the effect of the radion quantum 
fluctuations. 

To quantify the effect, we consider the power spectrum of the effective energy density, 



PE 



o ■ 



(5.1) 



Focusing on the positive tension brane, we define the power spectrum P n (k) of pe normalized by the background 
energy density p\ = 3H^_/k"^ as 



((Pe) 2 
(Pa) 2 



W) 



(5.2) 



where SVi-i is the surface area of the (n — l)-dimensional unit sphere; S n -\ — 2tt u ^ 2 /T(n/2). We find 

k n-lp /la 2TTK 2 n Hl- 2 (H+£) n - 2 , xn+3| ff (l) f-frnll 2 

k Fn[k) ~ 9(n~l)(n-2) 2 N {n _ 2)K coshz + { ^ ^("-sjM 



(H+l) 



71-2 



m— 3 



(-kr,)^o 9(n - l)(n - 2) 2 N (n _ 2)K cosh z+ r (n - 1/2) 2 ' 
Let us introduce the density parameter £Ie for the rms effective energy density of the radion field as 



PA 



(5.3) 



(5.4) 



which is obtained by integrating k n ~ 1 P n (k) over k and taking the square root of it. When £Ie is comparable 
to, or larger than the order of unity, the fluctuation of the brane is non-negligible. If we integrate the above 
over k, we encounter the divergence from large k. However, this ultraviolet divergence is the same as the usual 
one in the Minkowski background. Hence, to discard the divergence, we cut off the integration at the Hubble 
horizon scale; k = kn = The result is 



Sn- 



6(27t)"- 



-k n ~ 1 P{k) 



-i/\q\ «=4 135 y 7T cosh z+ (cosh z- — coshz+) 



(k 4 H + )(H + £). 



(5.5) 



Here and below we fix the dimensions to n = 4. Note that H + £ = sinh z + is the ratio of the AdSs curvature 
radius to the dS4 curvature radius, and k 2 = k 2 coshz + /£. If we introduce the 4-dimensional Planck scale in 
the single flat brane limit, £ v \ = ka/ ^/cosh z+, we have the well-known relation, l 2 ^ = £\j£- Below we assume 
£ v \ is the present-day Planck length in our universe. 

To discuss the cosmological const rain ts, it is convenient to introduce the 5-dimensional gravitational length 
scale £5 — k 2 / 3 . In terms of £5, Eq. (5.5) is re-expressed as 



135 V 7T y/cosh z- - cosh z + \£ 



3/2 



We see that the quantum radion effect will be negligible if 

coshz_ — coshz + > (H + £) 2 



3/2 



(5.6) 



(5.7) 
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Given the ratios H + £ and £s/£, this constrains the value of H- to be larger than a critical value, 



yA + (H-i) 2 > y/l + (H+£) 2 + (H + £f . (5.8) 

Since £5 is the fundamental scale of the 5-dimensional theory, it is natural to assume £5 < £; otherwise the 
whole system will be in the quantum gravitational regime. Note that (£ 5 /£) 3 = (£ P \/£) 2 - With this assumption, 
let us consider the following two limiting cases. 

(a) H+Kl: 



In this case, Eq (5.8) is satisfied for practically all values of H_ > H+ except for the range very close to 
that is, the constraint becomes 



(b) H+l > 1: 

In this case, we have 

f; - 1 £<" + '> (f ) 3/2 ■ 0"») 

Using the equality (£s/£) 3 / 2 = £ p i/£, the right hand side is rewritten as H + £ v \. For any reasonable scenario 
of inflation, we must have H + £ v \ <C 1. Hence the constraint on the value of H— is very mild in this case 
as well. 

Thus, with the assumptions that H + £ p \ <C 1 and {£^/£Y = (£ P \/£) 2 ^ 1, the constraints obtained in the above 
two cases can be concisely expressed as 

ln(§£) >(l+tf+£)^f. (5.11) 

To summarize, as far as the quantum radion fluctuations are concerned, the negative tension brane can be very 
close to the positive tension brane, even within a distance smaller than the AdSs curvature radius. 



VI. SUMMARY AND DISCUSSION 



In this paper, taking up the Randall-Sundrum type two-brane scenario in an (n + l)-dimensional spacetime, 
we have investigated the quantum fluctuations of the relative displacement of the branes, which we called 
the quantum radion. We have considered the cases of the flat two-brane and de Sitter two-brane systems 
simultaneously. Adopting the so-called Randall-Sundrum gauge, we have first solved the linear gravitational 
perturbation equations that describe the radion mode. 

Then assuming the perturbation of the form with a fixed z-dependence that solves the gravitational equation 
in the z-direction, where z is the extra-dimensional coordinate orthogonal to the branes, we have derived the 
effective action for the radion. With the effective action at hand, we have quantized the radion assuming the 
radion state is the Bunch-Davis vacuum. 

We have analyzed the effect of the quantum radion on the brane using the effective Einstein equations derived 
by Shiromizu, Maeda and Sasaki |y| in which the radion adds an effective energy-momentum tensor at the 
linear order in the field amplitude. Although the radion has the negative mass squared —nH 2 on the de Sitter 
brane where H is the Hubble parameter, we have found that the corresponding instability does not show up 
on the brane. We have noted, however, that the anisotropic stress induced by the radion is unusually large, 
though it still decays in time as a -1 for a fixed k, where a is the cosmic scale factor and k is the comoving 
wavenumber, and the rms value of the effective energy density for a fixed k decays as a~ 3 irrespective of the 
spacetime dimensions n. 

Focusing on the positive tension de Sitter brane with the Hubble rate H + , which models the braneworld 
inflation, we have estimated the rms total energy density of the quantum radion by integrating over k up to 
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aH + . We have introduced the density parameter f2 e which describes the relative magnitude of the radion energy 
density to the background energy density, and discussed the condition fls <C 1 on the model parameters. For 
H + £ p \ <§; 1 and £ 5 < £, which are reasonable to assume for the background to be not in the quantum regime, we 
have found that practically any choice of the location of the negative tension brane is allowed. This implies the 
quantum radion does not seriously affect the braneworld inflation scenario, at least at the linear perturbation 
order. 

The most intriguing question remained now is if the analysis here at the order of the linear perturbation 
level is sufficient. Naively, one may regard the effective action of the radion we have obtained as a piece to 
be added to the total effective action for an effective 4-dimensional theory that includes the gravity. Then the 
variation of the radion effective action with respect to the 4-metric will give the energy-momentum tensor of a 
scalar field with a negative mass squared, which would grow and diverge in contrast to the decaying effective 
energy-momentum tensor we analysed in this paper. However, such an energy-momentum tensor is quadratic 
in the radion, which is beyond the accuracy of our analysis at the linear-order perturbation level. We do need 
to investigate the higher order perturbation to truly see the effect of the quantum radion on the brane. This 
issue is left for future study. 
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APPENDIX A: EFFECTIVE ACTION OF RADION DERIVED AND COMPARED 



In this appendix, we give an outline of the derivation of the action given in Eq. (4.2) and compare i t wit h the 
effective actions of radion derived in [ ^o|Jl^1 . The action is obtained by substituting the results in Eq. ( |3.25|) into 
the action of the system given in Eq. ( |4.l| ). Note that 7^„ denotes the metric of the n-dimensional space-time 
with a constant curvature /C and the covariant derivative with respect to 7^. We do not demand ^(x^) to 



satisfy the field equation (3.26), so that the deg ree o f freedom associated with the brane fluctuation rema ins. 

We first substitute the results given in Eq. (3.25) into the bulk actions Ir_2A,i and -Tr-2A,2 of Eq. (4.1) 
Carrying out the calculation with the help of the commutation rule, 



(Al) 



results in vanishing of all the second-order terms exc ept for the terms that can be cast into the surface term on 
Yi±. Thus, to the second order in (p, the action (4.1) is found to have support only on Yl±: 




- V + \kfh + - hb- 2 U n ip - ( n - l)b- 2 lChv 



1 



-b 2 <j±k 2 (—ipn n ip — nJCip 2 



k^:=^b 2 d r (b- 2 h^). (A2) 



Note that ip here is defined as the displacement of the corresponding brane measured from the side of smaller 
values of r, that is, 



±L(z±;a)tp. 



(A3) 



Here, it is worth noting that the above action is different from the one obtained in || (see Eq. (A. 19)). 
This is because the above action 1^ is obtained by substituting a particular form of the metric to the action 
of the system, while t he a ction of || is constructed without any restriction on the metric. The particular form 
of the metric in Eq. (3.25) solves the equation of motion in z-direction, and thus already extremize the action 
in the z-direction. Therefore, is an action to be maximized in the x^-directions only, which is to give the 
equation of motion for hp V as a function of ip(x^). Meanwhile, the action obtained in H is an action to be 
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maximized in both the z— direction and ir^-directions independently to give the equations of motion for both 
and p Ipj . 
Let us define Ch and C v as 



C h := -k^h^ + kfh + h^ip^ - hb^Uip - (n - l)b- 2 tChp 
C v := ^h^p^ u - hb- 2 aip - (n - l)b- 2 ICh<pj + \^- 2 a ± n 2 



tpOtp — nJCp' 



(A4) 



Then is express as 



r(2) _ J_ 



We find that the variation of Ch with respect to gives the boundary condition of h^ v on the brane, while 
that of C v with respect to ip gives the equation of motion for tp. 

With the substitution of Eq. (3.25) to Ij, , Ch vanishes and only C v remains. Using the commutation rule 

d n x^( (&_cW 2) (z-) - b + d z u {2) (z+)) <p (□ + nK) 



(|Al|), C v is simplified to give 

r(2) _ r (2) ,(2) _ (n - 1)^ 



<P, 



which is Eq. (|4.2j). 

The effective action for the radion has been derived previously by Chiba p3] and by Chacko and Fox jl3| in 
the longitudinal gauge. However, in our opinion, there seems to exist some subtleties in their derivations that 
may be worth clarified, apart from the difference in the choice of gauge. We therefore compare their derivations 
with ours. 

First, we need to transform the metric perturbation in Eq. ( 3.25| ), which is given in the RS gauge, to that in 
the longitudinal gauge in which the locations of the branes are unperturbed, h rp = and /i M „ cx 7^. With the 
coordinate transformation of the form, 



r^r- L(z; a)p , x» -> x' 1 - {aU + ip {2) }p lf * , 
we obtain the metric perturbation in the desired gauge: 



-2{n-2)(b/£y- n <p, h 



0. 



V = 2b 2 (b/ey- n ^ v 



(A5) 
(A6) 



Let us first consider Chiba's effective action of radion given in Eq. (14) of |20[. The metric ansatz Chiba uses 
to derive the action is 



ds 2 = hldr 2 



r + f{x»)e 2kr . 



(A7) 



in which g^ is a general 4-dimensional metric that depends solely on x^. Here, the coordinate z in his paper is 
replaced with the coordinate r to avoid con fusion. Linearizing the metric (A7) with respect to /(x^), it reduces 
to the metric perturbation of the form ( |A6| ) with the following correspondence between his notation and ours: 



x» ^ ix^ . 



(A8) 



where we have fixed n = 4. (It may be noted that Eq. (A7) can be extended to the general n dimensions simply 
by replacing h = r + f(x fl )e 2kr with r + f {x^)e^ n ^ 2 ^ kr in the case of K. — 0.) With the above correspondence, 
we find his result, when linearized, does not agree with ours. 

Let us examine where the difference comes from by following his derivation. Substituting the decomposition 
of the Ricci scalar of the metric (A7), 



R[g] = -20k 2 + e 2kh R[g] + 6k 2 e 4kr+2kh h-\-l + 2kfe 2kr )V p fV p f + (total divergence term in x») , (A9) 
into his starting action of the system, Eq. (1) of EQ], and integrating it with respect to r, we obtain 
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e 



-3fcM 3 (a - l)(V (+) /) 2 - 3jfeM 3 (a~ 2 e~ 4(a_1)fc/ - 1) 



-<7( + ) / <f X^-g^) _Cr (-) / dx V~9(~) ' ( A10 ) 

where the following correspondence is understood: 

.9(±)^ <-> <7±^ > M 3 <->«T 2 , a^{i/b-f, cr (±) ^cr±, (A <-> A„+i/(2k 2 )) . (All) 



Comparing the action ( A10 ) with Eq. (5) in Q, it seems that the first term 5fcM 3 (a 2 e — 1) is 

missing there, which comes from the integration of the Ricci scalar of the background RS bulk spacetime, —20k 2 



in Eq. ([A9|). Furthermore, if the correspondence (All) we deduced is correct, the values of the brane tensions 



<7( ± ) = ±3M 3 fc as stated in Eq. (3) of Q are wrong by the factor 2; i.e., the correct values are a(±j = ±6kM 
Besides, in his starting action where the integral with respect to the extra dimension is given by twice the 
integral over r in the range (0, r c ), that is, from the positive tension brane to the negative tension brane, the 



geometrical boundary terms of the bulk spacetime, Ik,i and Ik. 2 as in our action of the system in Eq. (4.1) 
should be added. 



With these corrections we find 



S c (<T i±) -» ±6fcM 3 ) + - I K - 



= 2 / A,/=^{|L (i - le-^-w) e^R[g] - 3fcM 3 (a - l)(V (+) /) 2 } , 

where we have used the fact that Kau = —ke h g tiu . Since our effective action of the radion is obtained by 
fixing g M „ = ri^i, , the action to be compared is 



S c (a {±) -» ±6kM d ,g^ -> r,^) + I K .i - I K ,2 = 2 / (ficyf^i - 3fcM 3 (a - l)(V (+) /) 2 \ . (A12) 



Using Eqs. flAq ) and (All) and the correspondence N2K ^ {a — l)/2 which follows from them, we find that 
the above action agrees with our result. 

Let us next consider the derivation by Chacko and Fox and their result. The correspondence between their 
notation and ours are 

G M N^9ab, G^^q^u, M 3 «-> (4k 2 ) -1 , A b ^A 5 /k 2 , A «-> a+ , A^ <-> er_ , 

a^JJT 1 , ff«->l, f^b 2 , ^^2i 2 (p, g^^j^, (A13) 



from which we deduce 



dr'f- 1 <-> r^JVajc ■ (A14) 



Using the above correspondence, we find that their result given in Eq. (34) of |T^] exactly agrees with ours. 

There is, however, a slight subtlety which may be worth pointing out. It is the following. Chacko and Fox 
state that they work in the compact extra dimension with Z2 symmetry. Therefore there is no need to introduce 
the geometrical boundary terms in the action. In fact, the form of the linearized action given by Eq. (33) of (u| 
is clearly free from such terms. However, if we perform the integral over r of the linearized action, there appear 
contributions from the branes due to the curvature singularities. The interesting fact is that the contributions 
cancel out with those from the tension terms in the longitudinal gauge. Perhaps it is related with the fact that 
the coordinated locations of the branes are unperturbed in the longitudinal gauge. This implies that one may 
arrive at the correct answer by simply ignoring the singular contributions from the branes and simultaneously 
the tension terms, and then integrating the linearized action only over the bulk between branes. 

In the present paper, we took a different approach. That is, we divided the covering space of the bulk into 
the two patches bounded by the branes. This division requires the additional geometrical boundary terms for 
each patch to the action of the system. If we apply our method to their starting action, Eq. (1) of Jlj|, it is 
rewritten as 
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Scf = 2 / d*xV~G(2M*R - A B ) - I d A xy 



r=0 



+2 I d 4 xA 

lr—0 



lr=0 

lj (4M 3 K) - 2 



-GAo- / (PxV-GAa 



d 4 ^ 



r— a 
3 i 



-G (4M d iT) 



(A15) 



where AT^ is the extrinsic curvature of the r-constant hypersurface defined with the normal vector pointing 
toward the increasing direction of r. Apart from the difference in notation as expressed in the above, this 
form of the action is the o ne w e used as the starting point of our calculation. See Eq. ( f4.1| ). Substituting the 
perturbation given in Eq. (A6) into the above action, we find that the terms arising from the first integral in 
Eq. QA15D th at can be cast into the total divergence form cancel out with the perturbation of the last four terms 
of Eq. ( A15), i.e., the tension terms and the geometrical boundary terms. This cancellation gives the same form 
of action as the one in their Eq. (33) by restricting the integral only over the bulk, and consequently leads to 
their final result, Eq. (34). 



APPENDIX B: EFFECTIVE BRANS-DICKE PARAMETER 

In this appendix, we derive the effective Brans-Dicke parameter on the n-dimensional positive tension brane 
with curvature JC — 0, 1. 

The action for the Brans-Dicke theory in the n-dimensional space-time is 



I = 



1 



167rG r , 



qd n x 2A) - rVD^H* 



q d n xC 71 



(Bl) 



where q^ v is the metric of the n-dimensional space-time, wbd is the Brans-Dicke parameter, and C rn is the 
Lagradian density of matter fields. Let ^(x^ 1 ) = fo exp[M / (a; AI )]. Then, the equations of motion linearized in 
W(x») are 



+ Aq„ v = — - t m „ + [D^Du - q^Og] W{x»), 
Jo 



f a q w(x») 



n — 1 + (n — 2)w B d 



iirG n T - 2A/ ) , 



(B2) 
(B3) 



where t^ v is the energy-momentum tensor that arises from C m . Using Eq. (B3), Eq. (B2) can be put into the 
following form: 



n n — 1 + [n — 2)^bd 
_ 87rG„ (n - 2)w B D 



fa n — 1 + (n - 2)w B d 



d>irG n n — 1 

fo n—l + (n— 2)wbd 



D^D V - —q^U^ 



W(x») (B4) 



where t m „ := t^ v - (l/n)q^T. 

To obtain the effective Einstein equations on the positive tension brane for the energy density small compared 
with the tension of the brane, we substitute Eq. (3.8) of 0] into Eq. (2.2) of this paper. The result is 



w G +Aq = K *r + (w- 2 *% , V/b 



\n-2 r 



D„D 



u,^v Quis*— 'a 
n 



(B5) 



Comparing Eq. (B4) with Eq. (B5), we find that the effective gravity on the brane takes the Brans-Dicke 
form with the following identifications: 



87rG r . 



(n - 2)w B d 



k 8nG r , 



n — 1 



fo n-l + (n - 2)u B D ' 

Wb+r- 2 



£N, 



(n-2)K. 



2£N {n _ 2 )K fo n-l + (n- 2)w B d 



(B6) 
(B7) 



Eliminating fo from Eq. (B6), we obtain 
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^bd = (n - 1)- 



cosh ^/Kz + 



(sinh vXz + /vX)"- 2 



{n-2)K ) 



while eliminating wbd from Eq. (B£), we obtain 



(n — 2) cosh vXz + 



iV, 



(n-2)AC 



■ ( sinh VlCz + /VX 



For /C = 1 and n = 4, we have N/ n _ 2 )ic = cosh2:_ — coshz + and 



cosh; 



sinh 



-N. 



2K ■ 



8irG4 = 



K 2 t 



z+ 



2 cosh z + 



sinh 2 z + 



2K 



(B8) 



(B9) 



(BIO) 



For K — and n — 4, we have A r („_ 2 )/e = (z 2 — z 2 _)/2 and Eqs. (B8) and (B9) recover the result of Garriga 
and Tanaka B . 



APPENDIX C: E av A LA COSMOLOGICAL PERTURBATION THEORY 



In this appendix, we express the effect of the radion on the brane in the language of cosmological perturbation 
theory. We follow the notation of |21|. To begin with, we list a few important facts to be kept in mind: 



(a) The background is de Sitter space-time. 
3 ) T ^v = ~Ev/ K n 1S gauge-invariant as it is. 
(c) t^ v is traceless. 

The points (a) and (b) are somewhat related. As long as the background bulk spacetime is anti-de Sitter, 
the point (b) is always true. Independent of this fact, from the brane point of view, the matter perturbation 
variables are automatically gauge-invariant on the pure de Sitter space. For example, if we consider the density 
perturbation, the gauge transformation law for Sp is given by 



6p~^Sp = 5p+(n— l)aH(p + p)T , 



(CI) 



where T is the gauge function describing the shift in the time coordinate; 77 — > fj = 77 + T. Thus the density 
perturbation on exactly de Sitter space where p+p = is always gauge- invariant, and the same is true for all the 
matter perturbation variables. It should be noted, however, this special property of the de Sitter background can 
lead to an erroneous result if one tries to consider the perturbation in a gauge specified by conditions involving 
the matter variables. We shall come back to this issue later. For the moment, what we can claim is that the 
gauge-invariance is guaranteed provided one considers the perturbations in a class of gauges whose conditions 
involve only the geometrical variables such as the metric. Let us call this class of gauges the geometrical gauges. 
First we introduce the matter perturbation variables: 



<5T% = - P 5{ri) Y{x l ) , ST , = pq^Y^) , 



6T 



PttlM Y(x l ) + PTr T (ri)Y i j (x i ) , 



where Y(x l ), Yi(x l ) and Yij{x l ) are the spatial harmonics satisfying 
(V k V k + k 2 )Y(x l ) = 0, Y i (x i ):=-k- 1 V i Y{x i ), Y^x*) := AT 2 

and ST^v is just the effective radion energy momentum tensor, 

1 



.2 E V V ■ 



(C2) 



(C3) 



(C4) 



Note that the usual velocity perturbation, v = pq/(p + P), cannot be defined because of p + P = on the 
background. Also note that the traceless nature of t^ v implies hl = —5/(n— 1). As for the metric perturbation, 



we introduce the variables as 
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h^dx^dx" = -2a 2 AYdri 2 - 2a I BY J di 1 dx 3 + a 2 [2H L 8 t3 + 2H T Y lj ] dx'dx 1 , (C5) 

where we have chosen the spatially flat chart of the de Sitter space. 

As noted above, the choice of gauge is irrelevant as long as we choose a geometrical gauge. The perturbation 
equations are then derived by simply writing down the energy momentum conservation law ST^ = on the 
unperturbed background. Therefore, for definiteness, let us choose the Newton (or longitudinal) gauge in which 
the shear of the constant time hypersurfaces vanishes, i.e., k~~ 1 H' T — B = 0, and denote the matter variables in 
this gauge as 

6 = -(n-l)Tr L = A s , q = Q s , 7r T = n. (C6) 
Adopting the scale factor a as the time variable, we find 

d k 

a—A s +nA s = -Q s , (C7) 

da art 

d k ( 1 n — 2 \ 

a Q s+nQs= I As + n . (C8) 

da art \ n — 1 n — 1 / 

These equations can be combined to give a second order equation for A s , 

2 d 2 A . . d A / 1 k 2 , \ 7i-2 k 2 s 

We see from the above equation that if there were no anisotropic stress II, then A s cx a~ n or a~(" +1 ) after 
horizon-crossing. However, in the present case, we have II cx k 2 a~ 2 ip oc ka" 1 . This slow decay rate of the 



anisotropic stress acts as a source to the energy density A s . With II cx ka~ , the right hand side of Eq. (C9) 
behaves as k 3 a~ 5 , which implies A s cx fc 3 a~ 3 (for n > 3). The behavior of Q s is found as Q s oc k 2 a~ 2 . To be a 
bit more precise, the leading order behaviors of A s , Q s and II are given by 



C f k \ 3 „ C ( k \ 2 „ k 



A ' = - (n-l)(n-3) UJ ' °'~UJ ' "^W' (C10) 
wher e C is a constant. These results are of course fully consistent with those obtained in the main text, 



Eq. ( 4.15 ). Thus one may say that this unusual behavior of A s is caused by the large (though decaying) 
anisotropic stress. 

It is instructive to re-express the above result in the so-called comoving gauge in which the rj ^constant 
hypersurfaces are chosen in such a way that 

5T° j =pqY j = 0. (Cll) 

Note that this condition involves the matter variable q. A peculiarity of the comoving gauge is that the 
density perturbation in this gauge is not equal to A s , although the matter perturbation variables should be 
gauge-invariant on the pure de Sitter space-time as noted above. The density perturbation on the comoving 
hypersurface, which we denote by A, is related to A s as pl| 

A = A s + (n-l)^-Q s . (C12) 

One immediately sees from Eq. ( |C1(]| ) that the Q s term dominates and the leading order behavior of A is given 
by 

A = c ^r (C13) 

which decays much slower than A s . The cause of this seemingly inconsistent result is the comoving gauge 



condition (Cll) which forces q to be zero. As q is a gauge-invariant quantity, this condition would never be 
fulfilled and thus the comoving slice does not exist. In fact, it is straightforward to show that all the metric 
perturbati on va riables are ill-defined in the comoving gauge. Therefore, although the variable A is well-defined 



as in Eq. ( C12| ), this A does not represent a density perturbation in any gauge. 



Let us analyze the behaviors of the metric perturbation variables in the Newton gauge. Although not essential, 
for simplicity, we take Ht = which implies B — 0. We denote the metric variables in the Newton gauge by 
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A = * , H L =<I>, (C14) 
Then the (0, /Lt)-components and the traceless part of the (i, j)-components of the Einstein equations give 

a 

(n-3)$ + f = -^Pn = 3^j n. (C15) 
Inserting the leading behaviors of A and II given above to these equations, we find 

— < C16 > 

One may find it a bit surprising that these variables grow as a, seemingly indicating an instability. The resolution 
is that the Newton gauge is not really a good gauge. One should examine if there is a different choice of gauge 
in which the metric perturbations behave regularly, and there exists indeed such a gauge. By the shift of the 
time slice rj — > fj — i] + T, A and Hi transform as 

A^>A = A-H—(aT), H L -> H L = H L - HaT . (C17) 
da 

Then it is easy to see that the leading terms of $ and W are simultaneously eliminated by T = Q/(aH) =constant. 
Thus the above apparent instability is just a reflection of the bad choice of gauge. 
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